effects_ot = import(here("data","effects_ot.rds"))

fig_range = effects_ot |> 
  group_by(condition) |> 
  summarise(range = max(effect, na.rm = T) - min(effect, na.rm = T)) |> 
  arrange(desc(range)) |> 
  ggplot(aes(x = range)) +
  geom_density(fill = "grey", alpha = .5) +
  geom_vline(xintercept = 0.047, lty = 2, color = "red") +
  theme_classic() +
  annotate("label", x = 0.05, y = 10, label = "Average High/Low Income Difference (4.7%)", hjust = 0) +
  labs(x = "Range (max-min) of monthly CAMCE values over 2020 election cycle",
       y = "Density")

ggsave(here("figures","appendix","figure_s9.png"), fig_range,
       dpi = 600, units = "in", width = 6, height = 3)

ccat("See 'figure_s9.png' for Figure S9 output")